Spin dynamics in the 5 = 1/2 weakly-coupled chains antiferromagnet BaCu 2 Si207. 
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Inelastic neutron scattering is used to identify single-particle and continuum spin excitations in the 
quasi-one-dimensional S = 1/2 antiferromagnet BaCu2Si2C>7. In the data analysis, close attention is 
given to resolution effects. A gap in the continuum spectrum at A c = 4.8 meV is directly observed. 
The gap energy is in excellent agreement with existing theoretical predictions. Below the gap, the 
spectrum is dominated by sharp accoustic spin wave excitations. 
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I. INTRODUCTION 

Despite the apparent simplicity of the quantum S = 
1/2 antiferromagnetic (AF) Heisenberg model, the quest 
for a complete understanding of this system remains one 
of central problems in magnetism. Semi-classical spin 
wave theory (SWT) is known to provide an adequate de- 
scription of the 3-dimensional (3D) version of the model. 
After extensive theoretical and experimental studies it is 
fair to say that the purely one-dimensional (ID) case is 
rather well understood as well. Its static and dynamic 
properties totally defy the conventional spin wave pic- 
ture. In ID the ground state looks nothing like the 
classical two-sublattice Neel state, and is described as a 
"marginal spin liquid" , with power^jlaw spatial spin cor- 
relations and no long-range orders The spectrum con- 
tains no single-particle magnon excitation, but is instead 
a diffuse continuum described in terms of two-particle 
"spinon" statesflu Of great current interest is the quasi- 
1 dimensional case of weakly interacting quantum spin 
chains. Studies of such systems provide the missing link 
between the exotic quantum-mechanical behavior of the 
ID model and the more familiar semi-classical spin wave 
dynamics realized in three dimensions. 

A great deal of experimental information on the qua 
ID Heisenberg model was gained in studies of KCuFVl 
Particularly exciting was the recent observation of a novel 
longitudinal excitation, polarized along the direction of 
ordered moment, and thus totally absent in SWT.El The 
recently characterized BaCu2Si207 (Refs. |l0|-|i"2|) is an- 
other model weakly coupled chains compound, shown to 
be particularly useful for studying transverse spin fluctu- 
ations, i.e., those polarized perpendicular to the ordered 



moment. In a recent short paperS we have demonstrated 
that the transverse excitation spectrum in BaCu2Si207 
has a unique dual nature. At high energies the imaginary 
part of dynamic susceptibility is dominated by contin- 
uum excitations, characteristic of the quantum ID model. 
However, at energies below the characteristic strength of 
inter-chain interactions, the continuum is "cleaned up" 
and replaced by sharp single-particle states, similar to 
conventional spin waves. A further analysis suggested 
that continuum excitations appear only above a well- 
defined threshold energy. Unfortunately, for purely tech- 
nical reasons, this gap in the continuum could not be 
observed directly. In the present paper we report a de- 
tailed study of continuum excitations in the vicinity of 
the energy gap. 



II. AN EFFECTIVE SPIN HAMILTONIAN FOR 

BaCu 2 Si 2 7 

Before we proceed to describe the experimental proce- 
dures and results we shall briefly review the structure and 
relevant magnetic interactions in BaCu2Si207, that are, 
by now, rather well documented. The silicate crystallizes 
in an orthorhombic structure, space group Prima, with 
lattice constants arei-a = 8.862(2) A, b = 13.178(1) A 
and c = 6.897(1) A.LL3 Magnetic properties are due to 
Cu 2+ ions that form weakly-coupled antiferromagnetic 
chains running along the crystallographic c-axis. Cor- 
respondingly, the temperature dependence of bulk mag- 
netic .susceptibility follows the theoretical Bonner- Fisher 
curvecil with an in-chain coupling J — 24 meV.Ej A 
similar estimate was obtained in direct inelastic neutron 
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scattering measurements of the ID zone-boundary energy 
tkuzB of magnetic excitations, related to J through J 
2ftwzB/7r-l!^ This measurement gives J = 24.1 meV.tHl 
Nearest-neighbor spin-spin separation within the chains 
is equal to c/2, so the ID AF zone center qn — it cor- 
responds to I = 1, where Q — (h,k,l) denotes a wave 
vector in reciprocal space. 

Weak interactions between the chains result io-Jong- 
range antiferromagnetic ordering at Tn = 9.2 K.E2I The 
ordered moment is parallel to the chain axis and eX|teap^ 
olates to a saturation value of 0.15 /ib a t T — ► 0.E2TO 
In the ordered state the alignment of nearest-neighbor 
spins is parallel along the a, and antiparallel along the 
b and c directions, respectively. The 3D AF zone-center 
is located at Q = (0, 1, 1). The relevant inter-chain ex- 
change constants were determined from measurements of 
spin wave dispersion along several directions in reciprocal 
space.E3 The effective spin Hamiltonian was established: 

H = ^ ' JSij^ n Sij }7l +i 

i,j,n 

~t~ i.j ,n d" Jy & i,j i,j-\-l,n 

+ J 3 S (i) 

Here the sum is taken over all spins in the system. The 
indexes i and j enumerate the spin chains along the a and 
b axes, respectively, and n labels the spins within each 
chain. The experimental values of inter-chain exchange 
parameter are:t3 

J x = -0.460(7) meV, 
J y = 0.200(6) meV, 

2J 3 = 0.152(7) meV. (2) 

The spin chains in BaCu2Si2C>7 are not straight, but 
slightly zig-zag. The magnetic sites are displaced from 
high-symmetry (t, 0, j) positions as described in Ref. [l(]. 
For the neutron scattering measurements described be- 
low, this fact has one unfortunate consequence. The 
3D magnetic dynamic structure factor S(Q,u>) probed 
experimentally is not that of an ideal straight chain, 
Sq(Q,w), but related to it through the following rather 
complicated relation: 

S (Q, uj) = cos 2 (2tt/i(5 2; ) cos 2 (27rfc^) cos 2 (27rZ<5 z )S (Q, uj) 
+ cos 2 (2%hS x ) cos 2 (27rfc<5 y ) sin 2 (2tt 1 5 Z )S (Q + (100), uj) 
+ cos 2 (27TM,,) sin 2 (27rfc^) cos 2 (27r/<5 z )S' (Q + (011), uj) 
+ cos 2 (2ttM 2; ) sin 2 (27rfc^) sin 2 (27r/(5 z )5o(Q + (111), uj) 
+ sin 2 (2nhS x ) cos 2 (27rfc^) cos 2 (2 7 r/ ( 5 z )S' (Q + (001), w) 
+ sin 2 (27rM x ) cos 2 (27rfc^) sin 2 (27r/(5 z )5o(Q + (101), uj) 
+ sin 2 (27r/«y sin 2 (27rfc<5 y ) cos 2 (27r/(5 z )5o(Q + (010), uj) 
+ sm 2 (2nh5 x ) sm 2 (2nk5 y ) sin 2 (2tt I S Z )S {Q + (110), uj) (3) 

In this equation 5 X » 0.028, 5 y « 0.004 and 5 Z « 0.044 
are Cu 2+ displacements from high-symmetry positions 



as described in Ref. [10[|*| It is important to emphasize 
that Eq. || is an exact result. At small wave vectors 
only the first term is relevant and the dynamic structure 
factor is that of a straight chain. However, at higher wave 
vector transfers the contributions of other terms become 
significant. 

III. EXPERIMENTAL PROCEDURES 

The main goal of the present work was to reliably re- 
solve low-energy spin wave and continuum excitations 
near the ID AF zone-center in BaCu2Si2C>7. Previous 
studies indicated that in BaCu2Si2 07 the continuum sets 
in between 3.5 meV and 5.5 meV energy transferO The 
most interesting dynamic range is therefore between zero 
and 10 meV. The choice of spectrometer configurations 
in inelastic neutron scattering experiments was therefore 
governed by the conflicting needs of intensity, high wave 
vector resolution along the chain axis (required to over- 
come the steep dispersion along the chains) , energy reso- 
lution, and the need to cover a wide range of energy trans- 
fers. Since it appears impossible to simultaneously satisfy 
all these conditions, several different experimental setups 
with complementary characteristics were utilized. Quan- 
titative comparisons between data sets collected in differ- 
ent configurations was essential to verifying the validity 
of the results. In all cases we used the same BjaGu 2 Si207 
single-crystal sample as in previous studies .llillla In all 
configurations the crystal was mounted with the (0, fc, I) 
reciprocal-space plane in the scattering plane of the spec- 
trometer. 

Ideally, the measurements are done around Q = 
(0,0,1) (<7|i = 71"), where the extra terms in Eq. || are 
negligible, and all magnetic scattering is due to spin 
fluctuations transverse to the direction of ordered mo- 
ment (c axis). Unfortunately, at this wave vector, cer- 
tain geometric constrains restrict the accessible range 
of energy transfers to just over 4 meV. Measurements 
in this range were performed on the IN14 cold-neutron 
spectrometer installed at Institute Laue-Langevin (ILL), 
using a very high energy and wave vector resolution con- 
figuration (Setup I): guide — 40' — 40' — open collima- 
tions, Ef = 3 meV fixed final-energy neutrons, and a Be 
filter after the sample. To cover a wider energy-range 
and still have reasonable wave vector and energy reso- 
lutions we employed IN14 in "anti-W" geometry with 
guide — 40' — 60' — open collimations and final neutron 
energy fixed at Ef = 5 meV (Setup II). The data were 
collected around Q = (0,0,3) (q» — 3ir), where the con- 
tributions of both the first and second terms in Eq. are 



*We found that Eq. 8 in Ref. |lC] is incorrect. Equation g 
is the corrected version of this formula, rewritten to take S y 
into accounted. 
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significant. In this geometry measurements could be ex- 
tended up to 8 meV energy transfer. A similar configura- 
tion was used in combination with a horizontally-focusing 
analyzer (Setup III). While in this latter mode the wave 
vector resolution is considerably worse than with a flat 
analyzer, the intensity gain is quite substantial. Energy 
transfers of up to 17 meV were achieved using the IN22 
spectrometer at ILL, with guide — 60' — 60' — open colli- 
mations and final neutron energy fixed at Ef = 14.7 meV 
(Setup IV). In this case a Pyrolitic Graphite (PG) filter 
was used after the sample. 
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FIG. 1. Left: Typical constant-Q scans collected 
in BaCu2Si207 at the ID AF zone-center using different 
spectrometer configurations. Heavy solid lines represent a 
semi-global fit to the data as described in the text. Shaded 
areas are contributions of single-particle excitations given by 
Eqs. Dashed lines show the continuum contribution, as 
expressed by Eq. |l^. Arrows indicate the slight dip in the 
observed intensity that corresponds to the continuum energy 
gap A c . Right: Evolution of the calculated FWHM resolution 
ellipsoids in the course of the corresponding scans, plotted in 
projection onto the (I, %u>) plane. Solid lines represent the 
spin wave dispersion relation. 

In all experiments sample environment was a standard 
"ILL Orange" He-4 cryostat and all data were collected 
at either T = 1.5 K or T = 2 K. For the cold neutron 
experiments the background was measured in constant- 
Q scans, away from the ID AF zone-center, where no 
magnetic scattering is expected to occur, due to a step 
dispersion along the chain axis. In all cases the back- 



ground was found to be completely flat and featureless, 
typically 2.5, 4 and 5 counts per minute for setups I, 
II and III, respectively. The measured background was 
fit to a straight line and is subtracted from all the data 
sets shown below. For the thermal neutron experiment 
(Setup IV) the background was measured on an empty 
sample can and subtracted from the experimental data 
in a similar fashion. 
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FIG. 2. Constant-Q scans collected in BaCu2Si207 at the 
ID AF zone-center at different momentum transfers perpen- 
dicular to the chains using Setup III. Lines, shaded areas and 
graphical representations of the resolution functions on the 
left are as in Fig. ^. 



IV. RESULTS 

A. Transverse spin excitations 

Representative constant-Q scans measured at the ID 
AF zone-center qii = ir, 3ir using Setups I, II and III 
are shown in Fig. [I]. To the right of each scan we plot 
the calculated FWHM resolution ellipsoids in projection 
on the (Z, Tiuj) plane to show their size and orientation 
relative to the spin wave dispersion branches originating 
from the 1st and 2nd terms in Eq. (solid lines). Three 
distinct components of the excitation spectrum are ap- 
parent. The sharp peak at hu « 2 meV is the single- 
particle spin wave mode given by the first term in Eq. |^. 
This feature dominates the low-energy spectrum around 
Q = (0, 0, 1), where the other terms in Eq. |^ are negli- 
gible thanks to the small value of corresponding trigono- 
metric prefactors (Fig. [l] (a) and Figs. 1,2 in Ref. pd] ). 
At <7|| = 3tt however, the second term in Eq. || becomes 
significant and gives rise to the 3 meV peak in Figs. |l| 
(b) and (c). In essence, this second peak is an "echo" 
of the same spin wave excitation visible by virtue of a 
non-trivial 3D arrangement of magnetic sites. The third 
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feature in Fig. [I] is the diffuse scattering above the spin 
wave peaks that, as will be discussed in detail below, is 
not a resolution-enhanced "tail" of the spin wave peaks, 
but a separate entity. This scattering is attributed to 
continuum excitations in the system. In Figs. ^ (b) and 
(c) one can see a weak yet reproducible intensity dip that 
separates this continuum scattering from the underlying 
spin waves (arrows). 
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FIG. 3. Constant-Q scan collected in BaCu2Si207 at 
the ID AF zone-center using the thermal-neutron Setup IV. 
Lines, shaded areas and graphical representation of the reso- 
lution function on the left are as in Fig. [j]. 

Figure. || shows constant-Q scans collected with Setup 
III at different wave vector transfers perpendicular to the 
chains. The dispersion of the lower-energy spin wave is 
quite apparent. An almost total absence of dispersion in 
the 3 meV peak is consistent with that it originates from 
Q+(100) scattering (Eq. ||), and the known spin wave dis- 
persion along the (1, k, 1) reciprocal-space rod (Ref. [l^, 
Fig. 6). In all cases the continuum is observed, and the 
intensity dip occurs at roughly the same energy (arrows). 
There seems to be no obvious wave vector dependence of 
the continuum intensity along the b* direction. The en- 
ergy resolution of the thermal-neutron experiment (setup 
IV) is too coarse to resolve the two spin wave branches 
and the continuum. Figure. [| shows an energy scans col- 
lected in this configuration at Q = (0, 0, 3). 

Several constant-.^ scans measured with different con- 
figurations are shown in Fig. |]. At 3 meV energy trans- 
fer (Fig. |(a) ) two spin wave peaks are perfectly resolved 
using the highest-resolution configuration (Setup I). This 
proves that at these energies the excitation spectrum is 
entirely due to single-particle excitations. At the same 
energy transfer Setup II lacks the wave vector resolu- 
tion to resolve the spin wave peaks, as can be seen from 
Fig. |] (b). However, as will be discussed in detail in 
the next section, at higher energies (Fig. ^ (c,d)) while 
the resolution is sufficient to resolve individual spin wave 
peaks, only a single broad feature is observed. As ex- 
plained in Ref. [0], this behavior is a result of the exci- 
tation continuum "filling in" the space between the spin 
wave branches. The trend continues at even higher ener- 
gies, as can be seen from Fig. || (e) and (f), that shows 
scans collected using Setup IV. 



60 
40 
20 




(a) 


I Setup 1 




t » ftcu=3 meV 




11 Q=(0,0,/) 






2.9 3.0 3.1 3.2 3.3 
/(r.l.u) 

FIG. 4. A series of constant- .E scans along the spin chains 
in BaCu2Si207. Lines and shaded areas are as in Fig. [I]. 



B. Data analysis 

Resolution effects, particularly wave vector resolution 
along the chain axis, are clearly a major factor in the 
experiments discussed here. In order to reliably distin- 
guish between continuum scattering and the high-energy 
resolution "tails" of the single-particle peaks, one has to 
de-convolute the signal and the resolution function of the 
spectrometer. In practice this is done by constructing a 
parameterized model cross section, numerically convo- 
luting it with the calculated resolution, and adjusting 
the parameters to best-fit the data. Confidence in this 
procedure can be gained by simultaneously fitting data 



4 



sets collected in spectrometer configurations with differ- 
ent resolution characteristics, while using the same model 
cross section function and parameters. 



1. Model cross section function 

A simple yet accurate description of spin dynamics in 
weakly-interacting Heisenberg chains is provided by the 
Mean Field/ Random Phase Approximation (MF/RPA) 
treatment of inter-chain interactions.E30 Within this 
framework, the T = transverse-polarized spectrum of 
a weakly-ordered quasi-lD S = 1/2 system has two dis- 
tinct contributions. 
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FIG. 5. Measured intensity of transverse-polarized contin- 
uum scattering in BaCu2Si207, obtained by subtracting the 
simulated single-mode contribution from constant-Q scans as 
explained in the text. The solid line is a guide for the eye. 

The first contribution is from transverse-polarized 
single-particle excitations (spin waves), associated with 
the presence of long-range antiferromagnetic order. For 
BaCu2Si207 the spin .wave dispersion has been inves- 
tigated in some detai£3 and found to be in excellent 
agreement with MF/RPA theoretical predictions. In the 
present work we shall employ the same expression for the 
spin wave structure factor, as in Ref. [13: 
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In this formula Sgy L0 (Q,ui) and Sgy L0 (Q,u>) are single- 
mode (SM) dynamic structure factors for excitations 



polarized along the a and b axes, respectively, assum- 
ing straight spin chains in the system. The actual SM 
structure factors S^(Q,w) and Sg^(Q,uj) for the zig- 
zag chain in BaCu 2 Si20y are related to S§y L0 (Q,u}) 
and 5g^j (Q,cj) through Eq. ||[ The parameter A in 
Eq. represents an overall intensity -prefactor, while 
D x = 0.36 meV and D y — 0.21 meVlij are anisoropy 
gaps for the two spin wave branches. In Eq. [?] J(Q) 
is the Fourier transform of inter-chain coupling, and is 
given by 



J(Q) = 2J X cos(7r/i) + 2 J y cos(-7rfc) + 
2 J 3 cos(n(h + k)) + 2J 3 cos(n(h - k)). 



(8) 



The "mass gap" A in Eq. [7] defines the bandwidth of 
spin wave dispersion perpendicular to the chains. Its 
physical meaning deserves some comment. In the chain- 
MF model, in the magnetically ordered state, each spin 
chain experiences an effective staggered field gener- 
ated by the static staggered magnetization in adjacent 
chains. The staggered field produces a confining poten- 
tial between spinonsjn an S = 1/2 chain, yielding a two- 
spmon bound stateJl3 This new single-particle excitation 
has a gap A that scales as A oc iJ^ 3 .li3 When inter-chain 
interactions are switched on at the RPA level, the bound 
state acquires a dispersion perpendicular to the chain 
axis and takes the role of a conventional spin wave. At 
the 3D AF zone-center (the location of magnetic Bragg 
reflections) its energy goes to zero, as appropriate for a 
Goldstone mode. The mass gap A still enters the dis- 
persion relation, as in Eq. [?]. It is directly related to the 
strength of inter-chain interactions through:E3 



A ps 6.175IJ'! 



(9) 



From the known values of-inter-chain coupling constants 
for BaCu 2 Si 2 7 one getsM A = 2.51 meV. | 

The second contribution to the transverse dynamic 
structure factor in the MF/RPA model is what-remains 
of the 2-spinon continuum in isolated chainsOo Un- 
fortunately, a convenient analytical expression for this 
diffuse part of cross section, one that could be used to 



^It has recently been suggested— that the numerical coeffi- 
cient in Eq. [] may be incorrect Ji3 The actual value of this 
coefficient is actuaily irrelevant for this work, tndeed, the mi- 
croscopic exchange constants J1-J3 were deduced from the 
measured spin wave bandwidths in BaCu2Si2C>7 using the 
same vaiue of the coefficient. Note that it is A, and not 
the actuai J-values, that determine the spin wave dispersion 
perpendicular to the chains in Eq. |^. Unlike the exchange 
parameters, that characterize the model Hamiltonian, A is 
a true observabie physical quantity, equai to the measurabfe 
spin wave energy at Q — (0.5, 0.5, 1), where inter-chain inter- 
actions in BaCu2Si2C>7 cancef out at the RPA levef, and is 
thus independent of any theoreticaify-derived constants. 
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analyze the data, is not currently available. In our anal- 
ysis we instead employed a simple ansatz that has the 
same main characteristics as the MF/RPA result: i) At 
high energies the Muller-ansatz form that describes 2- 
spinon continuum scattering in isolated chainsE3 is re- 
covered; ii) The continuum has a sharp step-function 
low-energy cutoff (gap) A c = 2A at the ID AF zone- 
center q« = 7r; iii) Unlike the spin waves, the lower con- 
tinuum bound shows no dispersion perpendicular to the 
chain axis; iv) The gn- dependence of the lower bound 

is given by [A c ( 9 |,)] 2 = A 2 + £ J 2 sin 2 (g,,); and v) Ex- 
cept at "magic" reciprocal-space point (0.5, 0.5, 1), where 
inter-chain interactions cancel out at the RPA level, the 
continuum has no singularity at the lower bound. To an- 
alyze our data, where none of the measured scans were 
taken in the vicinity of the magic wave vector, we used 
the following empirical function to model the continuum 
part of the dynamic structure factor: 

s^yy(Q^) = aA [1 ~ c f 7Tl)] 

x 1 — 6[w- uj c (Q)} (10) 

V^ 2 -^ 2 sin 2 (<Z||) 

MQ)] 2 =A 2 + ^J 2 sin 2 (7r/) (11) 

The parameter a is the intensity of continuum scatter- 
ing relative to that of the spin waves, and 9(x) is the 
Heavyside step function. In essence, Eq. |l^ is simply the 
Miiller ansatz truncated at the expected lower continuum 
bound. For the 3D spin arrangement in BaCu2Si207 
the actual continuum cross sections S xx ' vv are related to 
S™ ' vv through Eq. |. 

2. Choice of adjustable parameters 

The goal was to perform a global analysis of all data 
sets collected in different configurations using as few ad- 
justable parameters as possible. All relevant exchange 
constants and anisotropy parameters have been measured 
previouslytj. These constants were fixed to the values 
quoted in Sec. [n| Since we have not performed absolute 
callibration of intensities measured in various setups, a 
separate value of the intensity prefactor A was used for 
each configuration. However, the same values of A were 
used for different scans measured using the same setup. 
In our analysis the continuum gap A c was used as a sep- 
arate variable, independent of A. As with the continuum 
intensity prefactor a, the same value of A c was used for 
all data sets measured in all configurations. 

In Eq. ||only the first two terms were found to be of any 
significance for the particular experiment, where the data 
were collected at small momentum transfers along the a- 
and b axes. Equation^ could thus be greatly simplified by 
setting the transverse zig-zag displacements S x and 5 y to 
zero. The relevant longitudinal displacement parameter 



S z can be expected to have a slight temperature depen- 
dence. Only the room-temperature value S z = 0.044 is 
known from structural data,t2l and in our analysis it had 
to be treated as an adjustable variable. The same value 
of S z was used for all data sets simultaneously. 

3. Fits to experimental data 

All together, 9 constant-Q and 8 constant-i? scans were 
analyzed, including all those shown in Figs. The 
model cross section, defined by Eqs. |7|-[Tl| and Eq. || was 
numerically folded with the calcualted Cooper-Nathans 
resolution functions^ for each setup. The calculation 
was based on the known geometry of the spectrometers 
and characteristics of the neutron guides, monochroma- 
tor and analyzer crystals, and collimators used in the 
experiments. The reliability of the resolution calculation 
was checked against the measured energy widths of elas- 
tic incoherent scattering from the sample in each setup. 
The folded cross section was scaled by the square of the 
magnetic form factor of Cu 2+ and fit to the data using a 
least-squares algorithm. 

The result of the fit (x 2 = 1-7) is shown in solid lines 
in Figs. [3]-||. The shaded areas are contributions of the 
single-mode part of the cross section. In the constant- 
energy scans the contributions of the spin wave branches, 
originating from the first two terms in Eq. 3d, are merged 
into a single profile. These two contributions are shown 
in shading of different intensity in the constant-Q data. 
The dashed lines represent the continuum part of the 
cross section. The following fit parameters were obtained: 
A c = 4.8(2) meV, a = 0.41(3) and S z = 0.041(1). 

It is gratifying to see that the model fit function re- 
produces the available data extremely well with just a 
few parameters. The refined value for S z is very close 
to its room-temperature value. The fitted continuum 
gap A c is within the error bar of the chain-RPA result: 
A c = 2A = 5.0 meV. 

To extract information on the excitation continuum 
without having to rely on a specific ansatz to describe its 
shape, we employed an alternative data analysis proce- 
dure. The data collected in constant-Q scans using Setup 
III below 3.5 meV were fit using only the SM part of the 
model cross section function (x 2 = 1.2). The result was 
then used to simulate the SMA contribution in the en- 
tire scan range and subtracted from the measured data. 
The remaining intensity is due to continuum scattering 
and is shown in Fig. pi where the results obtained at 
Q = (0,0,3), (0,0.5,3) and (0, 1,3) have been averaged. 
In this plot the threshold at around 2A meV becomes 
clearly apparent. 
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C. Attempt to observe the longitudinal mode. 



One of the most interesting predictions of the MF-RPA 
model for weakly-coupled S — 1/2 Heisenberg chains is 
the so-called longitudinal mode: a single-particle exci- 
tation polarized along the direction of ordered moment. 
This unique feature is totally absent in conventional spin 
wave theory, where the excitations represent the preces- 
sion of the order parameter around its equilibrium direc- 
tion, and are thus transverse-polarized. The longitudinal 
mode has been recently observed in the auasi-lD S = 1/2. 
system KCUF3 by means of unpolarizecB and polarizcdEj 
inelastic neutron scattering. In addition to the study of 
transverse-polarized excitations in BaCu2Si207, we per- 
formed a dedicated experiment to search for the longitu- 
dinal mode in this system. 
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FIG. 6. Reciprocal space/ energy- trajectory of the inelas- 
tic scan measured in BaCu2Si2 07 using setup V, showing the 
orientation of FWHM resolution ellipsoids relative to the spin 
wave dispersion. 

MF-RPA theory makes specific predictions for the en- 
ergy and intensity of the longitudinal mode. Unlike the 
transverse spin waves that have a substantial dispersion 
perpendicular to the chains and reach zero energy at the 
3D AF zone-center, the longitudinal mode is gapped and 
has very little dispersion in this direction— Its RPA dy- 
namic structure factor can be written as:l!3 

SSfoW.") = Al [1 - C0 2 S{nl)] 5y~ul(Q)], (12) 

MQ)] 2 - ^J 2 sin 2 H) 

+ ^[3|J'| + J(Q)], (13) 
7» 0.49 (14) 

From this equation we see that the longitudinal mode is 
expected to have a gap of \/3A and be approximately four 
times weaker in intensity than the transverse-polarized 
spin waves. 




ha (meV) 

FIG. 7. Constant-*^ scan collected in BaCu2Si207 at the 
ID AF zone-center and large momentum transfer perpendic- 
ular to the chains, using Setup V. The trajectory of the scan 
is as shown in Fig. |(| Solid line is an estimate of the con- 
tribution of spin fluctuation perpendicular to the direction of 
ordered moment. The dashed line shows the expected contri- 
butionr»f the longitudinal mode, as predicted by the MF/RPA 
model.kj The dotted line is the MF/RPA prediction for the 
total scattering intensity. 

To look for the longitudinal mode in BaCu2Si20y, 
measurements were performed at large momentum trans- 
fers perpendicular to the chain axis. In this geometry the 
angle between the ordered moment (c-axis) and the scat- 
tering vector is large, and both longitudinal and trans- 
verse excitations contribute to scattering. The data were 
collected at the NG5 cold 3-axis spectrometer at the Na- 
tional Institute of Standards and Technology Center for 
Neutron Research. The large momentum transfers per- 
pendicular to the chains enabled us to take full advantage 
of a horizontally-focusing PG analyzer. The chain axis 
was at all times aligned as closely parallel to the scat- 
tered beam as possible, ensuring that energy and c-axis 
wave-vector resolution widths remained at least as nar- 
row as in the standard flat-analyzer mode. The data 
were thus collected along a rather non-trivial trajectory 
in E — Q space, as illustrated in Fig. At all times the 
momentum transfer along the chains was maintained at 
the ID AF zone-center 1 = 1, Neutrons of Ef — 5 meV 
were used with a Be filter positioned after the sample and 
guide — 80' — 80' (radial) — open collimations (Setup V). 
The background (flat, 3.4 counts/ min) was measured by 
repeating the scan with the sample rotated by roughly 
20° around the a axis. 

The data are shown in Fig. |7[ At low energies, below 
3 meV energy transfer, we expect the main contribution 
to the scattering intensity to originate from the gapless 
transverse-polarized spin waves. After the inclusion of 
appropriate neutron polarization factors, the cross sec- 
tion for transverse excitations described above was fit to 
this low-energy portion of the scan (up to 3 meV trans- 
fer). All parameters, including A c and a were fixed at the 
values determined previously, and only the intensity pref- 
actor was treated as an adjustable variable. The result of 
the fit is shown in a solid line in Fig. We see that the 
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transverse-polarized cross section alone reproduces the 
measured scan rather well. The only inconsistency is the 
lack of a dip or even inflection point in the measured in- 
tensity at A c . However, the fit shows that experimental 
energy resolution should have been sufficient to resolve 
it. 



1DS=1/2 Heisenberg AF 



BaCu^O, 
Q=(0,1,/) 



BaCi^Sip, 
Q=(0,0,/) 




0.95 1.00 1.05 



•V* /(r.l.u.) 
FIG. 8. Re-distribution of spectral weight in a ID 5=1/2 
Heisenberg AF (a) that occurs when inter-chain interactions 
are switched on, as in BaCu2Si207 (b and c). 

In principle, one could subtract the fit from the ex- 
perimental data and claim that the residual intensity is 
due to longitudinal spin fluctuation. Given the severity 
of resolution effects, we were reluctant to push the data 
analysis to this level. What we did instead was simulate 
the expected contribution from the longitudinal mode, as 
given by Eq. |lj and an appropriate neutron polarization 
factor, and shown in a dashed line in Fig. [?]. Note that 
the "extra" intensity (arrow in Fig. [?]), not accounted for 
by the transverse dynamic structure factor, is centered 
at roughly the same energy as the expected longitudi- 
nal excitation. Moreover, the intensity "surplus", inte- 
grated in the range 3.5-5.5 meV (3.0(0.4) counts/(min 
meV)) is within the error bar of the expected intensity 
of the longitudinal mode integrated in this range (2.7 
counts/ (min meV)). If the longitudinal mode was a sharp 
single-particle excitation, as predicted by the RPA, a well 
defined peak centered at 4.2 meV should have been easily 
observed with our experimental statistics. This is illus- 
trated by the dotted line, which is the calculated sum of 
the transverse contribution obtained in the fit and that 
given by Eq. 
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for the longitudinal mode. The result can 
be reconciled with the RPA model, if one assumes that 
the longitudinal mode has a finite lifetime, i.e., an intrin- 
sic energy width of roughly 1 meV. This scenario would 
explain both the "extra" intensity and the absence of a 
sharp peak. It-is also consistent with experimental stud- 
ies of KCuF3,El where the observed longitudinal mode is 
broadened and has an intrinsic width of roughly 1/4 of 
its central energy. 



V. DISCUSSION AND CONCLUDING REMARKS 

Despite the obvious difficulties associated with wave 
vector- and energy resolution, the BaCu2Si2 07 data pre- 
sented above contain compelling evidence for a gap, or 
at least a pseudo-gap, in the transverse-polarized contin- 
uum. The gap value A c = 4.74.8 meV« 2A is in excel- 
lent agreement with predictions of the MF-RPA mpdel. 
Interestingly, non-linear spin wave theory (NSWT)t3 pre- 
dicts a much larger threshold energy. In this model the 
lowest-energy transverse-polarized continuum excitations 
are 3-magnon states, and a pseudogap occurs at roughly 
3A. An NSWT calculation for BaCu 2 Si 2 7 basecLan the 
known exchange parameters gives A c w 7.5 meV.cil 

Figure [s] is a cartoon illustration of how the spec- 
tral weight of transverse-polarized excitations gets re- 
distributed in a ID S = 1/2 Heisenberg AF when inter- 
chain interactions are switched on. The singularity on the 
lower bound in the ID system (Fig. ^ (a)) becomes sepa- 
rated from the bulk of the two-spinon continuum (Fig. || 
(b) and (c)). The spectral weight in its vicinity is con- 
solidated (arrows) to yield a sharp spin wave excitation 
(solid lines). 

Our results regarding longitudinal excitations are less 
conclusive. Clearly, additional measurements, possibly 
using polarized neutrons, will be required to fully resolve 
this problem for BaCu2Si207. The available unpolar- 
ized data shown above are, in fact, consistent with the 
MF-RPA model, but suggest * finite lifetime for the lon- 
gitudinal mode, as in KCuF3.1l 
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